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^ ■ ABSTRACT 
> ■ 

<^ : 

, We consider the effects of advection and radial gradients of pressure and radial drift velocity 

on the structure of accretion disks around black holes with proper description of optically 
I ■ thick/thin transitions. We concentrated our efforts on the models with large accretion rate. 

, Contrary to disk models neglecting advection, wc find that continuous solutions extending from 

' the outer disk regions to the inner edge exist for all accretion rates we have considered. We show 

. that the sonic point moves outward with increasing accretion rate, and that in the innermost 

disk region advection acts as a heating process that may even dominate over dissipative heating. 
Q ' Despite the importance of advection on it's structure, the disk remains geometrically thin. 

^ , Global solutions of advective accretion disks, which describe continuously the transition between 

^ I optically thick outer region and optically thin inner region are constructed and analyzed. 

• Subject headings: Accretion, accretion disks — black hole physics — hydrodynamics 

•i-H . 

^: 

' 1. Introduction 

The "standard accretion disk model" of Shakura (1972) and Shakura & Sunyaev (1973), that has been 
widely used to model accretion flows around black holes, is based on a number of simplifying assumptions. In 
particular, the flow is assumed to be geometrically thin and with a Keplerian angular velocity distribution. 
This assumption allows gradient terms in the differential equations describing the flow to be neglected, 
reducing them to a set of algebraic equations, and thereby fixes the angular momentum distribution of the 
flow. For low accretion rates, M, this assumption is generally considered to be reasonable. Since the end 
of the seventies, however, it has been realized that for high accretion rates, advection of energy with the 
flow can crucially modify the properties of the innermost parts of accretion disks around black holes. A 
deviation from a Keplerian rotation may result. 

Initial attempts to solve the more general disk problem only included advection of energy and the 
radial gradient of pressure in models with small values of the viscosity parameter, a = 10~^ (Paczyhski 
& Bisnovatyi-Kogan 1981), and it was shown that including radial velocity in the radial momentum 
equation would not change principally the results for such a small a (Muchotrzeb & Paczynski 1982). 
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Liang & Thomson (1980) emphasised the importance of the transonic nature of the radial drift velocity, 
and the influence of viscosity on the transonic accretion disk solutions was noted by Muchotrzcb (1983), 
who claimed that such solutions only existed for viscosity parameters smaller than a* ~ 0.02 — 0.05. 
Matsumoto et al. (1984), then showed that solutions with a > a* do in fact exist, but the nature of the 
singular point, where the radial velocity equals the sound velocity, is changing from a saddle to a nodal 
type and the position of this point is shifted substantially outwards in the disk. Matsumoto et al. (1984) 
also demonstrated the non-uniqueness of the solutions with a nodal type critical point for given Keplerian 
boundary conditions at the outer boundary of the disk (see also Muchotrzeb-Czerny 1986). Extensive 
investigation of accretion disk models with advcction for a wide range of the disk parameters. M and a, was 
conducted by Abramowicz et al. (1988), with special emphasis on low a. Misra & Meha (1996) considered 
optically thin two-temperature disk models and treated advection in the framework of the Keplerian disk 
model, but fixed the proton temperature somewhat arbitrarily at the outer boundary. Chakrabarti (1996) 
solved the advection problem containing shock waves near the innermost disk region, considering accretion 
through saddle points. Numerical solutions of accretion disks with advection have been obtained by Chen 
and Taam (1993) for the optically thick case with a = 0.1, and for the optically thin case (Narayan 1996). 
A simplified account of advection has recently been attempted, either treating it like an additional algebraic 
term assuming a constant radial gradient of entropy (Abramowicz et al. 1995; Chen et al. 1995; Chen 1995), 
or using the condition of self-similarity (Narayan & Yi 1994). 

Over the last few years it has become clear, that neglecting the advective heat transport at high M 
leads to qualitatively wrong conclusions about the topology of the family of solutions of the disk structure 
equations (see for example Abramowicz et al. 1995; Chen et al. 1995; Artemova et al. 1996). The disk 
structure equations without advection give rise two branches of solutions: optically thick and optically thin, 
which do not intersect if M < Mcr « (0.6 - 0.9)MEdd for a = 1 and Mbh = lO^Mg, where MEdd is the 
Eddington accretion rate (Artemova et al. 1996). For larger accretion rates there are no solutions of these 
equations extending continuously from large to small radii, and with Keplerian boundary conditions at the 
outer boundary of the disk (see also Liang & Wandel 1991; Wandel & Liang 1991; Luo & Liang 1994). It 
was argued by Artemova et al. (1996), that for accretion rates larger than Mcr advection becomes critically 
important and would allow solutions extending all the way to the inner disk edge also to exist for M > Mcr- 

The consistent transonic solutions for optically thick advective accretion disks were constructed in a 

wide range of values of the viscosity parameter a and mass accretion rate (sec Artemova ct al. (2001)) 
and character of singular points of obtained solutions was investigated. However for some values of mass 
accretion rate and viscosity parameter the corresponding value of optical depth in the inner region of 
accretion disk becomes less than unity and the used assumptions for optically thick disk approximation are 
violated. It requires to include into system of equations the transition formulae for radiative pressure and 
radiative flux (Artemova et al., 1996) to describe correctly the intermediate region between optically thick 
and optically thin zone. 

The goal of the present paper is to construct explicitly accretion disk models with optical depth 
transition for high M taking advective heat transport self-consistently into account. We also include radial 

gradients of pressure and radial drift velocity and we allow for the non-Keplcrian character of the circular 
velocity. Furthermore, we use the geometrically thin disk approximation because, the relative thickness of 
the disk in the considering models is less than unity. We show that solutions extending from large radii to 
the inner edge of the disk can be constructed even for accretion rates considerably larger than Mcr and 
the smooth transition between optically thick and optically thin regions in accretion disk appeares when 
the value of mass accretion rate becomes close to the critical one for considering models. We find that 
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advection is very important in the innermost disk region. 

In §2 we introduce our model and solution methods, while in §3 we discuss our results. 

2. The Model and the Method of Solution 

In this paper we will consider the full set of solutions to the disk structure equations with advection 
including the optically thick, optically thin and the intermediate zones in accretion disks. 



In our models we used two differents prescriptions for viscosity 

trci> = —aP , 

suggested by Shakura (1972) and the angular velocity gradient-dependent viscous stress, 



(1) 



(2) 



where u is the kinematic viscosity coefficient. 



We use from now on geometric units with G = 1, c = 1, use r as the radial coordinate scaled to rg — M, 
and scale all velocities to c. We work with the pseudo-Newtonian potential proposed by Paczynski and 
Wiita (1980), $ = —M/{r — 2), that provides an accurate, yet simple approximation to the Schwarzschild 
geometry. We normalize the accretion rate as m = M/MEdd, where M^dd = Lsdd = ^T^Mmp/uT, in our 
units. 

We use the same set of equations, ingredients and boundary conditions in our models as in the paper 
of Artemova et al. (2001)), except for changes required by the formulae for radiative flux and radiation 
pressure to describe correctly the intermediate zone in the accretion disk at high accretion rates. We used 
the same numerical method as Artemova et al. (2001) to solve this modified system of algebraic and 
difii'erential equations. 

The following equations are therefore modified: 

The vertically averaged energy conservation equation: 

Qadv ~ Q 1 
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are the advective energy, the viscous dissipation rate and the cooling rate per unit surface, respectively, T 
is the midplane temperature, k is the opacity, a is the radiation constant and tq is the Thompson optical 
depth. To = Kph. Here we have introduced the total optical depth to absorption, <IC tq. 
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and the effective optical depth 

T* = (roTa ) 

Where p is the density and h is the half-thickness of the disk. 

The equation of state for the matter consisted of a gas-radiation mixture is 

-Ptot = -Pgas + -Prad j (8) 

where the gas pressure is given by 

Pgas = pRT , 

where TZ is the gas constant. 

The expression for the radiation pressure is 

In our calculations we use the dimensionless parameter /3 to describe the relation of gas pressure to the 
total pressure in accretion disk 

p 



The specific energy of the mixture is 

pE = ^Pgas + 3Prad • (9) 

Our method allows us to construct a self-consistent solution to the system of equations from very large 
radii, r » 100, and down to the innermost regions of the disk. 



3. Results and Discussion 

Wc will now compare the solutions with advection and without it. In the "standard model", for 
accretion rates rh < rhcr = 36, (for a=0.5 and Mbh = IOMq, where rhcr is the critical accretion rate) there 
always exist solutions that extend continuously from large to small radii. When m > rhcr = 36 there are no 
solutions in a range of radii around r w 13, and therefore no continuous solutions extending from large radii 

to the innermost disk edge (see detailed discussion by Artemova et al. 1996, where however, the Newtonian 
potential was used, resulting in rhcr = 9.4 for q:=1.0 and Mbh = 10^ Mq). 

In Figure 1 and Figure 2 we plot the Thomson scattering depth tq and the effective optical depth 
r* as a function of radius, r, for different rh (for a=0.5 and Mbh = lOM©), clearly demonstrating that 
the solutions to the complete system of disk structure equations including advection and radial gradients 
have quite different properties at high m compared to the solutions of the standard disk model. For 
m < rhcr = 36 (dashed lines), including the gradient terms gives rather small corrections to the standard 
disk model. When m > 36 advection becomes essential and for rh > rhcr it changes the picture qualitatively. 
When rh > rhcr solutions do exist extending continuously from large radii to the innermost disk region 
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(the solid lines) where the solution passes through a "sonic point". However for small values of the mass 
accretion rate m < 0.1 all types of solutions (solutions without advection, optically thick solutions with 
advection and solutions with advection and optical depth transition) are very similar and have the same 
structure. 

Figure 3 shows the dependence of £in on accretion rate m [left panel], and the locations of the inner 
singular points (r^)™ as a function of m [right panel] for a = 0.5 and Mbh = lOM©. Dushed line 
corresponds to the advective optically thick solutions, and solid line - to the advective solutions with optical 
depth transition. 

Figure 4 shows the dependence of relation of gas pressure to the total pressure in accretion disk, Pg, 
on rh at the inner singular points for the same models as in Figure 3. The line styles are the same as in 
Figure 3. The change of value of P form 1 to corresponds to the change of a state from the gas pressure 
to radiative pressure dominated one. The thin disks with /3 ~ 1 are locally stable, whereas the parts of the 
disk in which /3 ~ could be thermally and viscously unstable at m ^ 100 (Pringle, Rees, & Pacholczyk 
1973). 

In Figures 5, 6, 7 we plot solutions with advection for m = 48 > rhcr, ot = 0.5, and Mbh = lOM©. We 
compare the solution obtained in optically thick approximation (dashed line), and the solution with optical 
depth transition (solid line). 

In Figure 5 the the smooth transition between optically thick outer region of accretion disk and 
optically thin inner region when the accretion rate is higher than the critical one is demonstrated (the solid 
line). 

Figure 6 shows the dependence of the temperature on the raduis for the same model as Figure 5 
and the transition from the outer regions with relatively low temperature to the inner edge of disk with 
temperature close to the l{fK is clearly seen (solid line). 

The corresponding dependence of density on the raduis is shown on Figure 7. When the optical depth 
transition formula is taken into account (solid line) the density of the disk decreases considerably in the 
inner regions of accretion disk, that corresponds to the optically thin high temperature state of the inner 
regions. 

Using analysis discussed by Artemova et al. (2001) we have determined a type of singular points in our 
numerical solutions. 

In the case of viscosity prescription (1) the solutions have two singular points. The inner singular 
points, {rs)in, can be saddles or nodes depending on values of a and m. But in the models discussed above, 
for a = 0.5 and m > 1.0, the inner singular point has always nodal-type. Note that the change of type from 
the saddle to the nodal one does not introduce any features in solutions. The outer singular points, (rs)out, 
are always of a saddle-type. In the case of viscosity prescription (2) the solutions have only inner singular 
points which are always of a saddle-type. 

We have obtained unique solutions for structure of advective accretion disk with optical depth 
transition in a wide range of accretion rates and a-parameters. Both viscosity prescriptions (1) and (2) 
have been investigated. The solutions corresponding to both prescriptions are very close for a ^ 0.1, and 
begin to differ at larger a. Due to the technical problems in calculation of the high viscosity models in the 
case of viscosity prescription (2) we consetrated our efforts on the study of high accretion rate models using 
prescription (1). 
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The main difference of the present study from the previous ones is in using more sophisticated numerical 
technique which accurately treats the regularity conditions in the inner singular point of equations. We have 
performed an analytical expansion at the singular point to calculate the derivatives of physical quantities. 
These derivatives help us to find the proper integral curve passing through the singular point. Theis 
approach allow us to avoid numerical instabilities and inaccuracies, appearing when only variables at the 
singular point, but not its derivatives, are included into a numerical scheme. 
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Fig. 1. — The dependence of the Thomson scattering depth on the radius for the models with a = 0.5 and 
Mbh = WMq. Dashed hnes correspond to the solutions without advection and m < rhcr = 36. Dotted 
lines correspond to the non-physical solutions without advection for rh = m,cr = 36 and m = 50 (from the 
center to the edge of the picture respectively). Solid lines correspond to the solutions with advection and 
the mass accretion rate higher than the critical one. Thick solid line corresponds to m = 36.0, and the thin 
solid line to m = 50.0. 
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Fig. 2. — The dependence of the effective optical depth on the radius for the same models as in Fig.l. The 
styles of the curves are the same as in Fig.l. 
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Fig. 3. — The specific angular momentum £in (left panel) and the position of the inner singular points (right 
panel) as a function of the mass accretion rate M for viscosity parameters a = 0.5 and Mbh = lOM©. 
The solid lines represent models with the optical depth transition and advection, whereas the dashed lines 
correspond to advective optically thick ones. 
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Fig. 4. — Ratio of the gas pressure to the total pressure at the inner singular point, (3s, as a function of the 
mass accretion rate M. Models with a = 0.5 and Mbh = lOM© are shown. The styles of the curves are the 
same as in Fig. 3. 
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Fig. 5. — The dependence of the effective optical depth on the radius for Mbh = WMq, rh = 48, a = 0.5. 
Dashed Une corresponds to the solution with advection and without optical depth transition, solid line 
corresponds to the solution with advection and optical depth transition. 
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Fig. 7. — The dependence of the density on the radius for the same models as in Fig.5. The styles of the 
curves are the same as in Fig.5. 



